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Dynamics and constraints of the dissipative Liouville cosmology 
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In this article we investigate the properties of the FLRW flat cosmological models in which the 
cosmic expansion of the Universe is affected by a dilaton dark energy (Liouville scenario). In 
particular, we perform a detailed study of these models in the light of the latest cosmological data, 
which serves to illustrate the phenomenological viability of the new dark energy paradigm as a serious 
alternative to the traditional scalar field approaches. By performing a joint likelihood analysis of 
the recent supernovae type la data (SNIa), the differential ages of passively evolving galaxies, and 
the Baryonic Acoustic Oscillations (BAOs) traced by the Sloan Digital Sky Survey (SDSS), we put 
tight constraints on the main cosmological parameters. Furthermore, we study the linear matter 
fluctuation field of the above Liouville cosmological models. In this framework, we compare the 
observed growth rate of clustering measured with those predicted by the current Liouville models. 
Performing a \ 2 statistical test we show that the Liouville cosmological model provides growth rates 
that match sufficiently well with the observed growth rate. To further test the viability of the models 
under study, we use the Press-Schechter formalism to derive their expected redshift distribution of 
cluster-size halos that will be provided by future X-ray and Sunyaev-Zeldovich cluster surveys. 
We find that the Hubble flow differences between the Liouville and the ACDM models provide a 
significantly different halo redshift distribution, suggesting that the models can be observationally 
distinguished. 

PACS numbers: 98.80.-k, 95.35.+d, 95.36,+x 
Keywords: Cosmology; dark matter; dark energy 



I. INTRODUCTION 

The comprehensive study carried out in recent years 
by the cosmologists has converged towards a cosmic ex- 
pansion history that involves a spatially flat geometry 
and a recent accelerating expansion of the Universe (see 
[lHgl and references therein). From a theoretical point 
of view, an easy way to explain this expansion is to con- 
sider an additional energy component, usually called dark 
energy (hereafter DE) with negative pressure, that domi- 
nates the Universe at late times. The simplest DE candi- 
date corresponds to the cosmological constant (see |9l-tll| 
for reviews). Indeed the so-called spatially flat concor- 
dance ACDM model, which includes cold dark matter 
(CDM) and a cosmological constant (A), fits accurately 
the current observational data and thus it is an excellent 
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candidate to be the model which describes the observed 
Universe. 

Nevertheless, the identification of A with the quan- 
tum vacuum has brought another problem: the estimate 
of theoretical physicists that the vacuum energy density 
should be 120 orders of magnitude bigger than the mea- 
sured A value. This is the "old" cosmological constant 
problem (CCP) Q. The "new" problem [l(j] asks why is 
the vacuum density so similar to the matter density just 
now? Many solutions to both problems have been pro- 
posed in the literature J12H14J . An easy way to overpass 
the above theoretical problems is to replace the constant 
vacuum energy with a DE that evolves with time. Nowa- 
days, the physics of DE is considered one of the most 
fundamental and challenging problems on the interface 
uniting Astronomy, Cosmology and Particle Physics and 
indeed in the last decade there have been theoretical de- 
bates among the cosmologists regarding the nature of this 
exotic component. 

In the original scalar field models [15| and later in the 
quintessence context, one can ad- hoc introduce an adjust- 



ing or tracker scalar field (f> [16[ (different from the usual 
SM Higgs field), rolling down the potential energy V{<p), 
which could resemble the DE [lfl[H ElHiij • This class 
of DE models have been widely used in the literature 
due to their simplicity. Notice that DE models with a 
canonical kinetic term have a dark energy EoS parameter 
— 1 < w^ < —1/3. Models with (w^ < —1), sometimes 
called phantom DE [2l|, are endowed with a very exotic 
nature, like a scalar field with negative kinetic energy. 
However, it was realized that the idea of a scalar field 
rolling down some suitable potential does not really solve 
the problem because <j> has to be some high energy field of 
a Grand Unified Theory (GUT), and this leads to an un- 
naturally small value of its mass, which is beyond all con- 
ceivable standards in Particle Physics. As an example, 
utilizing the simplest form for the potential of the scalar 
field, V(4>) — m^(f) 2 /2, the present value of the associ- 
ated vacuum energy density is pa = (V (</>)) ~ 10 -11 eU 4 , 
so for (</>} of order of a typical GUT scale near the Planck 
mass, Mp ~ 10 19 GeV, the corresponding mass of is 
expected in the ballpark of m^ ~ iJ ~ 10 -33 eV. 

Notice that the presence of such a tiny mass scale in 
scalar field models of DE is generally expected also on the 
basis of structure formation arguments |22h24J ; namely 
from the fact that the DE perturbations seem to play 
an insignificant role in structure formation for scales well 
below the sound horizon. The main reason for this homo- 
geneity of the DE is the flatness of the potential, which 
is necessary to produce a cosmic acceleration. Being the 
mass associated to the scalar field fluctuation propor- 
tional to the second derivative of the potential itself, it 
follows that 7710 will be very small and one expects that 
the magnitude of DE fluctuations induced by <j) should be 
appreciable only on length scales of the order of the hori- 
zon. Thus, equating the spatial scale of these fluctuations 
to the Compton wavelength of <p (hence to the inverse of 
its mass) it follows once more that m^ < H o ~ 10~ 33 eV. 

Despite the above difficulties there is a class of viable 
models of quintessence usually called dilatonic that are 
based on supersymmetry, supergravity and string-theory 
and which can protect, for specific potentials, the light 
mass of quintessence (for a review see Ref. |25| and ref- 
erences therein). In particular, in string theory, gauge 
and gravitational couplings are connected with the vac- 
uum expectation value of scalar field called dilaton (j> [26( . 
In this context, it has been proposed by some of us a 
specific model for the DE. This DE model being associ- 
ated with a rolling dilaton field that is a remnant of this 
non-equilibrium phase described by a generic non-critical 
(Liouville) string theory [27h3C| . We call this scenario ei- 
ther Q-cosmology or Liouville cosmology. This pattern 
is based on non-critical (Liouville) strings [3ll - l33l | , which 
offer a mathematically consistent way of incorporating 
time-dependent backgrounds in string theory. We note 
in passing that the presence of time dependent dilaton 
fields at late eras of the Universe, that characterizes Q- 
cosmology, may lead to phenomenologically interesting 
extensions of minimal supergravity models. The latter 



predict less dark-matter relic abundances than the con- 
ventional ACDM model, thereby allowing for more room 
for supersymmetry in collider (such as LHC) tests of the 
models (30,113. 

In this article we investigate the main dynamical prop- 
erties of the Q-cosmological model from the point of 
view of current astrophysical constraints, extending non- 
trivially earlier analyses [35, 36]. Initially, a joint sta- 
tistical analysis, involving the latest observational data 
[SNIa 3, H{z) [13 and BAO UlSj] is implemented. 
Secondly, we attempt to compare the matter fluctua- 
tion field of the Q-cosmology with that of the ACDM 
models by computing the growth rate of clustering. Fi- 
nally, by using the Press-Schechter formalism [41| and re- 
cently derived functional forms for the dark matter halo 
mass function [42], we show that the evolution of the 
cluster-size halo abundances is a potential discriminator 
between the Q-cosmology and ACDM models. We would 
like to stress here that the abundance of collapsed struc- 
tures, as a function of mass and redshift, can be accessed 
through observations [4J]. Indeed, the mass function of 
galaxy clusters has been measured based on X-ray sur- 
veys [44M46| , via weak and strong lensing studies |47H49| , 
using optical surveys, like the SDSS l50Ll5ll . as well as, 
through Sunayev-Zeldovich (SZ) effect [52j. In the last 
decade many authors have found that the abundance of 
collapsed structures is affected by the presence of a dark 
energy component |53l465l |. 

The structure of our paper is as follows. The basic el- 
ements of the Liouville cosmological model are presented 
in Section UH where we also introduce the cosmologi- 
cal equations for a flat Friedmann-Lemaitre-Robertson- 
Walkcr (FLRW) geometry. In Section ILUl a joint statis- 
tical analysis based on SNIa, H(z) and BAO is used to 
constraint the Liouville model free parameter. The issue 
related with the effective DE equation of state parame- 
ter and thus with the cosmic expansion is presented in 
Section IIVI The linear growth factor of matter pertur- 
bations is discussed in Section [V] while in Section I VII 
we discuss and compare the corresponding theoretical 
predictions regarding the evolution of the cluster abun- 
dances. Finally, the main conclusions are summarized in 
Section EH] 



II. THEORETICAL BACKGROUND ON 
Q-COSMOLOGY 

In non-critical Q-cosmological models the (effective, 
four-dimensional) Friedman equation is modified, includ- 
ing, apart from the standard matter and dark energy 
contributions denoted collectively by p, also Liouville- 
string off-shell corrections, Ap, which are not positive 
definite in general (without affecting, however, the over- 
all positive-energy conditions of this non-equilibrium the- 
ory): 
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The detailed form of the system of dynamical equa- 
tions, one of which is ([T]) is given by Ref. [66| in the 
Einstein frame [31j . i.e. in the frame where the scalar 
curvature in the (off-shell) target space effective action 
assumes the canonical Einstein form to leading order in 
the Regge-slope a' expansion, are given by 
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and q is the deceleration q = —da /a 2 . Above, an over- 
dot denotes derivative with respect to the FLRW cosmic 
time, t, p m , p m denote "matter" energy and pressure 
densities respectively, whilst 
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denote the dilaton energy density and pressure. In (Q|, 
the potential has the form V a ii = 2Q 2 exp (2</>) + V , 
where Q = Q{t) is the so-called central charge deficit 
function of the non-critical theory, which plays the role 
of a "gravitational friction" coefficient. The term V rep- 
resents higher-string loop corrections in the target space- 
string effective action, which are in general not known in 
a closed form. In our work the precise form of V will not 
play a role, at least for late eras of the Universe, where 
such quantum string corrections are sub-dominant. The 
variation of the central charge deficit Q(t) with the cos- 
mic time t is provided by the following consistency equa- 
tion, to leading order in the Regge-slope (a') expansion: 
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dt a 

The reader should notice that, although we have assumed 
a (spatially) flat Universe, the terms on the r.h.s. of ©, 
which manifest departure from criticality, act in a sense 
like curvature terms as being non-zero at certain epochs. 
The continuity equation for matter, which follows by 
combining the set of equations @, reads: 
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To proceed further one needs to make some extra assump- 
tions. First we assume that the matter-energy density is 
split as 



Qd + Qr + QS ■ 



(7) 



The first term on the right-hand-side of ([7]) refers to as 
"dust", Wd = 0, and includes the baryonic matter and 
any other sort of matter, characterized by w = 0, which 
does not feel the effect of the non-critical terms. Follow- 
ing Ref. [66| we may assume that dust does not couple 
to the dilaton and the non-critical-stringy corrections at 
late eras, so that its respective density, ga, obeys the 
standard energy conservation equation (c.f. ((6|)): 
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dt 
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The second term in ||7J) refers to radiation w r = | and 
the third term to an unknown sort of exotic matter which 
is characterized by an equation of state (hereafter EoS) 
having a weight w§. In our analysis for late epochs of 
the Universe, such as the eras of structure formation, we 
may assume that this sort of matter which would oth- 
erwise feel the non-critical terms in the evolution of the 
Universe, feels predominantly only the time-dependent 
dilaton <j> terms in ([6]), i.e. its density, g$, obeys approx- 
imately the following energy conservation equation: 
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dt 
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for late eras of the Universe, e.g. corresponding to red- 
shifts z < 2 , where analytic treatments of the Liouville 
Q-cosmology are available. In general the exotic matter 
equation of state w$(z) depends on the redshift, but in 
our analysis in this work we consider it to be a constant. 
For late epochs of the Universe, an approximate scaling 
behavior of the dilaton field may be found by making the 
further assumption that its configuration is not affected 
much by the presence of matter. In such a case, one may 
use the purely gravitational solution of Ref. |2q|. which 
expresses the dilaton in the form 1 : 



0~ ci + <f> ln(a(t)) 



(10) 



in units of the present scale factor o,q = 1, where ^>n, c\ 
are constants. In the purely gravitational theory of [28J, 
0o = —1, but here we want to be more general, in order 
to account for the presence of matter, and so we shall 
consider the constant ^>o as a phenomenological param- 
eter to be determined by fitting the data. However, we 



1 A similar behavior of the dilaton is also valid in the early Universe 
as discussed in |34|. Such forms depend crucially on the form 
of the dilaton potential, which should be of exponential form 
V ~ e~ cr ^, c some constant. Dilaton Q-cosmologies are indeed 
characterized by such potentials [351 I66H . 



shall assume 4>q < 0, which assumes a perturbatively 
weak string coupling g s — e* as the cosmic time elapses, 
which is crucial for the validity of our analytical treat- 
ment at late epochs of the Universe. In particular, this 
implies that string loop corrections, although non negli- 
gible, and in fact their inclusion may even lead to some 
interesting physical effects (see below) (35|, nevertheless 
they can be treated perturbatively, and thus can lead to 
analytic solutions at late eras. 

With the above in mind, we thus observe that from 
([ID]) and © we obtain: 



dgs 
~dt 



■Ml [ 1-y +w s (l + <l>o))gs^O 



(11) 



implying the following scaling behavior of the exotic mat- 
ter with the scale factor: 
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In particular, it has been found that 5 ~ 4 [36j can fit 
the BAO data [38]. Such a value is in the ballpark of 
theoretically expected values. For instance, the simplest 
phenomenological assumption on the pertinent EoS pa- 
rameter ws — — of the exotic dark matter fluid, consis- 
QS ' 

tent with Big-Bang-nucleosynthesis constraints, is that 
ws ~ 0.4 66]. To obtain S ~ 4, while maintaining the 
feature 4>$ < 0, one must use either values of \<fio\ <C 3, 
implying a weak dependence of the dilaton on the scale 
factor, or <f>o ~ — 1 and more generally the condition 



energy has been proposed by Sahni et al. [67[ • We would 
like to stress once again that the above formulas are valid 
for late eras, such as the ones pertinent to the supernova 
and other data (0 < z < 2). 

In this framework, the various fli contain contribu- 
tions from both dark energy and matter energy densi- 
ties. Specifically, H3 does not merely represent ordinary 
matter effects, but also receives contributions from the 
dilaton dark energy. In fact, the sign of O3 depends on 
details of the underlying theory, and it could even be neg- 
ative. In a similar vein, the exotic contributions scaling 
as (1 + z) 5 are affected by the off-shell Liouville terms 
of Q-cosmology. It is because of the similar scaling be- 
haviors of dark matter and dilaton dark energy that we 
reverted to the notation fjj, i = 2,3,(5 in Refs. [3a. l36j. 
In order to distinguish ordinary matter from dilaton dark 
energy effects that scale similarly with the redshift, one 
would have to perform also measurements of the effective 
dark energy EoS parameter (see Section IIVI) . 



III. LIKELIHOOD ANALYSIS 

Let us now discuss the statistical treatment of the ob- 
servational data used to constrain the Liouville model 
presented in the previous section. To begin with, we 
consider the Union supernovae la set of Kowalski et al. 
[4j which contains 307 entries. The likelihood estimator 
is determined by a XsNia statistics: 



in which case ©, (fTTj) become independent of w$. 

On the other hand, the BBN constraint is also satisfied 
in models where the dilaton dominates the early eras of 
the Universe [34]], in the absence of non-critical string 
effects, with the exotic matter behaving as dust ws = 
coupled to the dilaton. In this case, the scaling exponent 
(TT2"]) becomes: 5 ~ 3 + \(f>o\. The phenomenologically 
likely [36j value S ~ 4 is then achieved for \c/>o\ ~ 1, 
which is in the range of the analytical solutions for late 
eras found in [281 ] . and for which, as already mentioned, 
the energy equation (|11[) is largely independent of ws ■ 

Following the above discussion, we therefore 
parametrize the normalized Hubble parameter of 
this model, E(z) — H(z)/Hq, where z is the redshift, as 
follows [351 ]: 

E(z) = [q 3 (i + zf + n s (i + zf + n a (i + zf] 1/2 (13) 

with 
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It is interesting to mention that the current normalized 
Hubble function Eq. flTBl has a form which is close to a 
polynomial one. Note that a quadratic polynomial dark 
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where Zi is the the observed redshift, /i is the distance 
modulus fi = m — M = blogd^ + 25 and d^ is the lumi- 
nosity distance, 



dh{z) = (1 + z)r(z) r(z) = c 



dy 
H(y) 



(16) 



In addition to the SNIa data, we also consider the 
BAO scale produced in the last scattering surface by the 
competition between the pressure of the coupled baryon- 
photon fluid and gravity. The resulting acoustic waves 
leave (in the course of the evolution) an overdensity sig- 
nature at certain length scales of the matter distribu- 
tion. Evidence of this excess was recentl y fo und in the 
clustering properties of SDSS galaxies [38l - l40l ] and it pro- 
vides a suitable "standard ruler" for constraining dark 
energy models. We would like to remind the reader that 
in the Liouville cosmology the ordinary matter compo- 
nent which appears in the normalized Hubble parameter 
Eq. (fTBf is mixed with the "exotic" contributions. Thus 
we can not use the measurements of BAO derived by 
[38 - l40l . To overcome this problem Mavromatos & Mit- 
sou [361 ] have proposed a different estimator which is valid 



for Liouville cosmologies. This is 
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where zbao = 0.35 
B = 0.334 ± 0.021 
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Note that the measured value is 
Therefore, the corresponding 
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Finally, a very interesting geometrical probe of dark 
energy is provided by the measures of H(z) [33] from the 
differential ages of passively evolving galaxies (hereafter 
H{z) data). Note that the sample contains 11 entries 
spanning a redshift range of < z < 2. In this case the 
corresponding Xh function can be written as: 
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where H{z) is the Hubble parameter 2 , H(z) = HqE(z). 
In order to put tighter constraints on the correspond- 
ing parameter space of our cosmological model, the above 
probes are combined through a joint likelihood analy- 
sis 3 , given by the product of the individual likelihoods 
according to: Aot = £sNia x £bao x £h, which trans- 
lates in the joint \ 2 function in an addition: Xtot = 
XsNia + ^bao + Xh- The resulting best fit parameters for 
different values of o, are presented in Table Q] The current 
statistical results are in very good agreement with those 
found by |36| . 



TABLE I: Results of the likelihood function analysis. The 1 st 
column indicates the Liouville model. 2 n column presents 
the values of S used. 3 rd , 4 th and 5 th columns show the best 
fit parameters and the reduced Xtot- I n the final column one 
can find various line types of the models appearing in Figs, [l] 
and [3 
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-1.43 ±0.14 


0.27 ±0.02 
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solid 


(h 


4.1 


-1.64 ±0.16 


0.40 ± 0.04 
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dashed 
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-1.93 ±0.19 


0.59 ±0.06 


1.027 


dotted 
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-2.45 ±0.25 


0.95 ±0.10 


1.025 


triangles 


Qs 
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-3.29 ±0.33 


1.63 ±0.17 


1.024 


open cirles 



The reader may worry that negative dust-like energy 
densities indicate an instability of the model and an ob- 
vious violation of the positive energy conditions, and 
thus its immediate exclusion. However, as discussed in 
[35j, |36[ , and already mentioned in the previous section, 
the "dust" -like contributions, ^3, do not merely repre- 
sent ordinary matter effects, but also receive contribu- 
tions from the dilaton dark energy. In fact, the sign of 
f2 3 depends on details of the underlying theory, and it 
could even be negative. For instance, as argued in |35| . 
string loop corrections could lead to a negative Q3. In 
addition, Kaluza-Klcin graviton modes in certain branc- 
inspired models [69( also yield negative dust contribu- 
tions. In a similar vein, the exotic contributions scaling 
as (1 ± z) 5 are affected by the off-shell Liouville terms of 
Q-cosmology. 



IV. THE RECENT EXPANSION HISTORY 

In Fig. [T]we plot the Hubble parameter of the current 
Liouville models as a function of redshift, which for red- 
shifts z > 1.5 appears to be different in amplitude with 
respect to the corresponding ACDM model expectations 
(dot-dashed line). 
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FIG. 1: Comparison of the observed (solid circles [371 ]) and 
theoretical evolution of the Hubble parameter, H(z), using 
Ho = 73.8 ±2.4 km/s/Mpc HI]. The different Liouville mod- 
els are represented by different symbols and line types (see 
Table [J for definitions). For comparison, the dot-dashed line 
corresponds to the traditional ACDM model, while the thick- 
line error bars correspond to la //o-uncertainties. We do not 
plot the la i/o-uncertainties for the Liouville models in order 
to avoid confusion. 



The Hubble constant is H = 73.8 ± 2.4 km/s/Mpc 16811 . 
Likelihoods are normalized to their maximum values. In the 
present analysis we always report la uncertainties on the fitted 
parameters. Note also that the total number of data points used 
here is iVtot = 319, while the associated degrees of freedom are: 
dof= 317. Note that we sample ^3 G [-4, 1] and Cl s £ [0.1, 2] in 
steps of 0.001. 



As we have already mentioned in Section |TT] the density 
related with the f^ parameter varies as Qz{z) ex (l±z) 3 . 
Owing to the fact that the density of ordinary matter 
(baryons and usual cold dark matter), denoted g m from 
now on, follows the same power-law as ^3(2), Q m {z) oc 
(1 ± z) 3 , we are completely justified to assume that the 
usual matter density is one of the components of the den- 
sity Qzi. z )- In other words the Q m (z) — g m (z)/g c (z) den- 
sity parameter (where g c (z) = 3H 2 (z)/8wG is the critical 



density parameter) is included in the ^(z) parameter. 
On the other hand, as a result of the existence of non- 
trivial dilaton couplings to some dark matter species, as 
mentioned previously, there is dark matter with exotic 
scaling £ls(l + z) s , with 5 given by (fT2j) . Although this 
is not a dark energy contribution, nevertheless it is the 
dilaton couplings that modify its scaling, and therefore 
in this sense, the origin of this exotic scaling may be at- 
tributed to the existence of the interaction of dilaton dark 
energy with dark matter. However, since the fl$ > rep- 
resents massive matter, it may cluster and lead to non 
trivial contributions to structure growth in the Universe. 
In this sense, it is distinguishable from dark energy con- 
tributions that are non clustering. One of our points in 
this work is an attempt to distinguish phenomenologi- 
cally whether the current astrophysical data point out 
towards clustering of this exotic matter, and hence con- 
tribution to structure formation and growth, or not. 

To this end, we first write the normalized Hubble pa- 
rameter Eq. (Tr"5|) as follows 
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with 



ah(z) 2 = n«(i + z) s + (Q 3 - n m )(i + zf + n 2 {\ + zf. 

(21) 
Notice here that in the above formula we have been care- 
ful to isolate the matter contributions, pertaining to dust 
fid (for baryons and dust dark matter f2, = gio/ g c o), from 
the exotic matter terms, tt$ with non-trivial equation of 
state ws, corresponding to a scaling exponent S (fT2"|) . It 
should be stressed that the last term AH 2 (z) of the nor- 
malized Hubble function (|20[) encode the Q-cosmology 
corrections to the standard FLRW expression. 

Since the corresponding components exhibit smooth 
scaling with the redshift, the absorption of these extra 
effects to an "effective dark energy" contribution, with a 
non-trivial EoS is possible. In general, it is well known 
that one can express the effective dark energy EoS pa- 
rameter in terms of the normalized Hubble parameter 
701 
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3 din a 



1 - Q m (a) 



where in our case 



fl m (a) 



E 2 (a) 



(22) 



(23) 



Differentiating the latter and utilizing Eq. (|2"2")l we find 
that 



dQ,„, 

dln< 



3wuE(a)ft m {a)[l-n m (a)} . (24) 



After some simple algebra, it is also readily seen that 
the effective dark energy EoS parameter is given by (see 



w DE {a) 



wde(z) = - 



1- 
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din AH 2 
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din a 


L + 


1 
3 


din AH 2 


d\n(l + z) 



(25) 



(26) 



where a = (1 + z)~ 1 . From the above analysis it becomes 
clear that any modifications to the EoS parameter are in- 
cluded in the second term of Eq. ([26)) . A similar analysis 
for other DE models can be found in [7l], UM ■ 

In our case, on inserting Eq. (|2"Tj) into Eq. (121)1) it is 
straightforward to obtain a simple analytical expression 
for the effective dark energy EoS parameter: 



wde(-z) = — 1+ 



i 3(fi 3 - n m )(i + z) + 5Q s (i + z Y~ 2 + 2 ^2 



3 (n 3 - fim)(i + z) + n«(i + zy 



+ n 2 • 

(27) 



Notice that we will use the above quantity in Section fVl 



It thus follows that in the current cosmological context, 
the Liouville scenario as proposed by Ellis et al. [35| can 
be treated by an additional effective fluid with EoS pa- 
rameter defined by Eq. (|27|) . In order to investigate the 
behavior of the effective EoS parameter we have to know 
a priori the value of f2 m . The matter density fl m remains 
the most weakly constrained cosmological parameter. In 
principle, fl m is constrained by the maximum likelihood 
fit to the WMAP [6] and SNIa Q data in the context 
of the concordance ACDM cosmology. However, in the 
spirit of this work, we want to use measures which are 
completely independent of the dark energy component. 
An estimate of f2 m without conventional priors is not an 
easy task in observational cosmology. Nevertheless, var- 
ious authors, using mainly large scale structure studies, 
have attempted to set constraints on the f2 m parameter. 
In a rather old paper, Plionis et al. 73], using the motion 
of the Local Group with respect to the cosmic microwave 
background, found that O m ~ 0.30. From the analysis of 
the power spectrum, Sanchez et al. 74| obtained a value 
U m ~ 0.24. Moreover, the authors of Refs. [73, UM ana_ 
lyzed the peculiar velocity field in the local Universe and 
obtained the values J7 m ~ 0.30 and ~ 0.22 respectively. 
In addition, the authors of Rcf. [73], based on the clus- 
ter mass-to-light ratio, claim that Vi m lies in the interval 
0.15 — 0.26 (see also (78(] for a review). Therefore, there 
are strong independent indications that Q m lies in the 
range 0.2 < Q, m < 0.3. In order to compare our results 
with those of the flat ACDM model, we shall restrict our 
present analysis to an indicative value of fl m = 0.28. 

In order to visualize the redshift dependence of the ef- 
fective EoS parameter, we compare in the upper panel 
of Fig. [2] various flat Liouville cosmological models (see 
Table [l| . One can divide the evolution of the cosmic ex- 
pansion history in different phases on the basis of the 
varying behavior of the Liouville and ACDM models. We 



will investigate such variations in terms of the decelera- 
tion parameter, q(z) = — [1 — dhiH/dhi(l + z)], which is 
plotted in the lower panel of Fig. [2j We can divide the 
cosmic expansion history in the following phases: 

• At early enough times z ~ 2 the deceleration pa- 
rameters of both models are positive with qg < q^, 
which means that the ACDM model is more decel- 
erating than the Liouville model. 

• For 1.6 < z < 1.9 the deceleration parameters are 
both positive with qQ > q\, which means that the 
cosmic expansion in the Liouville model is more 
rapidly "decelerating" than in the ACDM case. 
Note that the effective EoS parameter of the Q- 
model tends to its maximum value, wde(z) ~ 1.4, 
while we always have u>de = —1 for the ACDM 
model (dot-dashed line); 

• Between 1.1 < z < 1.6 the deceleration parameters 
remain positive (with qQ < qa) but qQ decreases 
rapidly. 

• For 0.7 < z < 1.1 the traditional ACDM model 
remains in the decelerated regime (q\ > 0) but 
the Liouville model is starting to accelerate (<?q < 
0, wde < — l/3). At z ~ 0.9 the corresponding 
effective EoS parameter cross, for the first time, 
the phantom divide (wde = — 1) and it stays there 
for some time (0.2 < z < 0.9). 

• Prior to the present epoch the effective EoS pa- 
rameter crosses the phantom divide for the second 
time (z ~ 0.2) and then it remains close to wde — 
—0.90. The deceleration parameters are both neg- 
ative and since qQ > q\, the ACDM model pro- 
vides a stronger acceleration than in the Liouville 
model. At the present time we find qQ : o — —0.33 
and q\,a — —0.60 



V. THE GROWTH FACTOR AND THE RATE 
OF CLUSTERING 

In this section, we discuss the basic equation which 
governs the behavior of the matter perturbations on sub- 
horizon scales and within the framework of any dark en- 
ergy model, including those of modified gravity ("geo- 
metrical dark energy" ) . 

The basic assumption underlying the approach is that 
whatever we call "dark energy " contribution in our effec- 
tive approach, pertains to a substance that is non clus- 
tering to participate significantly in the growth of cosmic 
structures. In our Liouville cosmology, this issue is sub- 
tle. As we have discussed above, we do have two kinds of 
matter in the model, with two different equations of state 
and scaling, an ordinary dust, and an exotic dark matter 
component, with non trivial scaling exponent <5 (|12p . as 
a result of its coupling to the dilaton. The phenomeno- 
logically likely value S ~ 4 makes the behavior of exotic 
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FIG. 2: Expansion history. In the upper panel we display the 
evolution of the dark energy effective EoS parameter. The 
different Liouville models are represented by different sym- 
bols and line types (see Table [J for definitions). In the lower 
panel we compare the deceleration parameters of the Q2 and 
the concordance ACDM (dot-dashed line) models. The decel- 
eration parameter is plotted only for the Q2 model in order 
to avoid confusion. 



dark matter resembling (to a good approximation) that 
of relativistic matter (radiation), which does not cluster. 
We shall come back to this point later one. For the mo- 
ment let us restrict ourselves to the ordinary matter as 
the dominant kind that contributes to structure forma- 
tion via its clustering properties. 

For these cases, a full analytical description can be in- 
troduced by considering an extended Poisson equation 
together with the Euler and continuity equations. Con- 
sequently, we consider here the evolution equation of the 
matter fluctuations for models where the DE fluid has a 
vanishing anisotropic stress and it is not coup led to or- 
dinary matter (noninteracting DE see |72l |79h84| ) . It is 
well known that for small scales (smaller than the hori- 
zon) the dark energy component (or "geometrical" dark 
energy) is expected to be smooth and thus it is fair to 
consider perturbations only on the usual matter com- 
ponent of the cosmic fluid [85|- This assumption leads 
to the traditional equation for matter perturbations (see 
Appendix [5} : 



6 m + 2HS m = 4nGp m Sr, 



(28) 



where p m is the matter density which clusters in order to 
form cosmic structures. 

We would like to emphasize here that Eg. (|2"8l can be 
used also in the case of a dilaton field, as it has been 



8 



shown by Boisseau et al. [86| 4 In other words, Eq. (|28|) 
is valid also in the Liouville cosmology, as far as ordinary 
dust matter (including potential dark matter components 
scaling like dust) is concerned. 

Within this latter framework, we assume that for small 
scales (smaller than the horizon) the effective dark en- 
ergy component AH(z) appeared in the modified Fricd- 
mann equation (|20j) is smooth (homogeneous) which im- 
plies that Eq. ([25)1 is valid also in the Liouville model. 
If we change the variables from t to In a (-£ = H ,f ) 
then the time evolution of the mass density contrast (see 
Eq. ([28])) takes the following form js3| 

(ln<5 m )" + (ln<5 m )' 2 + (\n5 m ) X(a) = ~fi m (a), (29) 

where 

X(a) = ~-~w DE (a)[l-n m (a)]. (30) 

The prime denotes derivatives with respect to In a. Now, 
for any type of dark energy an efficient parametrization 
of the matter perturbations is based on the growth rate 
of clustering [88] 



/(«) = 



d In dm 
din a 



n^(o) 



(31) 



7 is the so called growth index [7l|, l79l I82J 1871 l89j| which 
plays a key role in cosmological studies, especially in 



4 The authors of 86] have worked in the so-called Jordan frame, in 
which the dilaton couples to the Einstein sacalar curvature term 
in the (four space-time dimensional) action through a function 

F(4>): 



S- 



I 



^m 



R- 



g^d^d^-VW + C^ 



st matter 



(ff/w) 



In this formalism there is an effective gravitational "constant", 
G e ff((j>), which is a function of the dilaton that replaces G in 
(128 D - Notice that the matter action does not couple to the 
dilaton in the Jordan frame. In our Liouville string cosmol- 
ogy case we have an additional component to the matter ac- 
tion, £ cxo ti c matter (flu^i $)> involving the exotic dark matter, 
with equation of state w& , and its interactions with the ordinary 
baryonic or dark matter with dust scaling. This depends on the 
dilaton in the Jordan frame. In the Liouville Cosmology we work 
in the Einstein frame, as already mentioned, by having redefined 
appropriately the gravitational field, in such a way that, in terms 
of the new metric field, there is a canonically normalized Ein- 
stein term in the effective action. Such a redefinition leaves the 
perturbative scattering amplitudes of low-energy string theory 
invariant, and hence is an allowed transformation in the string 
framework that does not affect the physical predictions. Un- 
der this redefinition, the matter action involves coupling of the 
dilaton with the dust-like matter fields, as a consequence of the 
redefined gravition. Nevertheless, we still assume for our pur- 
poses that such a coupling is very weak, and the only dominant 
couplings of the dilaton to matter fields in the Einstein frame are 
the ones pertaining to the exotic dark matter components. As 
we shall discuss below, such assumptions can be put directly to 
experimental tests by studying galactic growth data. 



light of recent large redshift surveys, such as the WiggleZ 
|90T r9lj . Indeed, it has been proposed that measurement 
of the growth index can provide an efficient way to dis- 
criminate between modified gravity models and DE mod- 
els which adhere to general relativity. As an example, it 
was theoretically shown that for DE models within the 
general relativity framework, the growth index 7 is well 
fitted by 7gr ~ 6/11 [82, 89]. Notice, that in the case 
of the braneworld model of Ref. [95[ we have 7 rs 11/16 
(see also [82J). while for the f(R) gravity models we have 
7(2) w 0.41 at z = H3. 

Now, inserting Eq. (13"TT) into Eq. (|29p and using simul- 
taneously Eq. (j2~4")l we obtain (see also [871 ]) 

3iUDE(a)O m (a)[l-fi m (a)]— h/ 2 + /X(a) = -O m (a) 

n (32) 
In this framework, Nesseris & Perivolaropoulos [89] used 

for the first time growth-rate data (most of the data are 

included in our Table ILT|) in order to constraint 7. They 

found (see also [82|, |87|) that a good approximation of 

the growth index, especially at z ~ 1, is given by: 



7 



3(w DE - 1) 
6wde - 5 



(33) 



Of course the above formula boils down to 6/11 for the 
usual ACDM [wde(z) = —1] model as it should. 5 

Since the effective EoS parameter of the Liouville cos- 
mological model varies strongly with redshift (see Fig. [2]) 
one would expect that the corresponding growth index 
would be a function of redshift. Following the notations 
of Polarski & Gannouji [96| we consider a first order ex- 
pansion of 7 in redshift 



li z ) = 7o +liz, 



(34) 



where 70 is given by Eq. (l33|) . Using Eq. (j34|) in (|3"Tj) and 
(|32]l we obtain (see also [96]) 



7i 



WX + 3(70 



!)WDE, (1 - Sim) - f^r 70 + \ 



In fL 



(35) 

where 70 is given by Eq. (|33[) . We would like to remind 
the reader that in the current article we use iude,o — 
—0.90 and ft m ~ 0.28 which imply that (70,71) — 
(0.550, -0.048). We also find that the latter set of values 
remains quite robust within the range of the physically 
accepted values of 51 m . Indeed as an example in the case 
of n m = 0.30 we find (70,71) ^ (0.549, -0.046). 



5 We should point out at this stage that in order to compare the 
observed growth history of the Universe with predictions coming 
from modified gravity models, which however is not our case, one 
cannot use the observed / t, s (z) but rather a combination of the 
growth rate of structure and the rms fluctuations of the linear 
density field on scales of 8h~ 1 Mpc, namely f(z)ag(z). This is a 
model independent quantity (see e.g. |9ll and references therein). 




the best fitted values for the parameters (O3 and fig) 
obtained in Section [TT] for each Liouville model. We then 
compute the corresponding consistency between models 
and data (Xmin/H) an d place these results in Table Hill 
For the Liouville models we find Xmin/H — 1-19 while 
in the case of the concordance ACDM cosmology we find 

Xmi n /H - °- 60 - Note that in the last column of Table Hill 
we list the corresponding results by excluding from the 
statistical analysis the observed growth rate of clustering 
at z = 3 [98(]. If we consider Q m — 0.30, we then find 
either xLJ 11 - I- 4 ? ° r X 2 min (z < 3)/10 ~ 0.70. 



FIG. 3: Comparison of the observed, [solid circles; see Ta- 
ble [II], and theoretical evolution of the growth rate of clus- 
tering f(z). Note that the different Liouville models are rep- 
resented by different symbols and line types (see Table U for 
definitions). The dashed line corresponds to the ACDM cos- 
mology. 



In Table UH we quote the observational values of the 
growth rate of clustering, with the corresponding error 
bars. This is an enriched version of the dataset used 
in 18 9f. which includes in addition data from the Wig- 
gleZ [9(| . The observed growth rate of clustering is given 
by fobs — fib, where (5{z) is the redshift space distor- 
tion parameter and b{z) is the linear bias. Observation- 
ally, one can estimate the /3(z) parameter by using the 
anisotropy of the correlation function. The linear bias 
factor can be defined as the ratio of the variances of the 
tracer (galaxies, QSOs, etc) and underlying mass density 
fields, smoothed out at the scale 8ft. -1 Mpc, at which 
the variance is of order unity: b(z) = as,tr(z)/as(z), 
where crgjr {z)Js measured directly from the sample. As 
discussed in [89j, the weak point of the fobs(z) data is 
the fact that the as(z) is defined at each redshift us- 
ing a particular reference (fiducial) ACDM model (see 
third column in Table ILT1) . Moreover, the growth mea- 
surements are correlated with the shape of the power 
spectrum, which makes the growth rate sample rather 
heterogeneous [9lL l92j . In our analysis, we ignore the 
effects of such a heterogeneity. Such simplifying assump- 
tions in treating the growth rate data have been used 
extensively in the literature in order to put constraints 
on the growth index 7, especially in the context of the 
braneworld cosmological model [93]. We would like to 
stress that in the present work, due to the above caveats, 
we do not use the fobs(z) data to constrain the cosmolog- 
ical parameters or the growth index. Instead, by using 
the constraints found from the cosmic expansion data 
(see section HTI) . we shall only compare the evolution of 
the predicted Liouville growth rate of structure with the 
data (for a similar analysis see [94|). 

In Fig. [3]we display the predicted growth rate Eq. (f3l]) 
together with the observed f b s (z) [97| . We compare the 
growth rate of clustering between data and models via 
a x 2 statistical test. For the model predictions we use 



TABLE II: Data of the growth rate of clustering. The cor- 
respondence of the columns is as follows: redshift, observed 
growth rate, cosmological parameters use by different authors 
and references. 



z 


Jobs, Re f 


(Clm,Ref,Cr8,Ref) 


Refs. 


0.15 


0.49 ±0.14 


(0.30, 0.90) 


[99 101] 


0.22 


0.60 ±0.10 


(0.27,0.80) 


[90] 


0.35 


0.70 ±0.18 


(0.24,0.76) 


[102] 


0.41 


0.70 ±0.07 


(0.27,0.80) 


[90] 


0.55 


0.75 ±0.18 


(0.30, 1.00) 


[103] 


0.60 


0.73 ±0.07 


(0.27,0.80) 


[90] 


0.77 


0.91 ±0.36 


(0.27,0.78) 


[99] 


0.78 


0.70 ± 0.08 


(0.27,0.80) 


[90] 


1.40 


0.90 ±0.24 


(0.25,0.84) 


[104] 


2.42 


0.74 ±0.24 


(0.26,0.93) 


[105] 


3.00 


1.46 ±0.29 


(0.30,0.85) 


[98] 



TABLE III: Statistical quantification of the comparison be- 
tween the observed growth rate of clustering (see Table [TTJ) 
and the theoretical model expectations. The last column cor- 
responds to comparison excluding the observed growth rate 
of clustering at z — 3 [98J . 



i/ll Xmin(^ < 3)/10 



Q-Data 1.19 
A-Data 0.60 



0.90 
0.45 



Finally, let us now examine the case in which the den- 
sity gs{z) oc (1 + z) s participates in the clustering. The 
situation is subtle. As we have discussed above, the ex- 
otic matter obeys a modified energy equation, as a re- 
sult of its direct coupling to the dilaton sources in the 
conservation Eqs. (j6]) and ((9]). Thus, in addition to the 
existence of pressure in the exotic dark matter fluid, due 
to a possibly non-trivial equation of state, w$, one has 
the effects of the dilaton coupling that also contribute to 
the growth equation. For 0o ~ — 1 in (J10[) . or more gen- 
erally <j> ~ —H, the modified conservation equation (fTTj) 
becomes independent of the equation of state Wg: 



gs + 3(1 ± v)Hq s ~ 0, 



(36) 
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where v = 1/3. This energy equation resembles that of a 
radiation fluid, with effective scaling law for the density 
qs ~ a~ 4 . However, the reader must bear in mind that 
the energy equation (|36l) stems from the time-dependent 
dilaton-sources modified energy equation ([9J, which for 
the case at hand, <f> ~ —H, is independent of the true 
equation of state ws . Within the standard Cosmological 
framework, the growth equation for radiation is modified, 
as compared to that for matter (|28|). by term s larger by 
a factor of about (1 + t/)(l + 3i/) = 8/3 [Toel ]: 



(37) 



,l+3v 



at 



6 S + 2HS S ~ 4ttG(1 + i/)(l + 3v)p s 8 s 

In this case the growing mode scales as 5$ ~ 
large redshifts which implies that the growth rate takes 
a large value /,5^1 + 3^ = 2. The latter value seems to 
be ruled out by the observed growth rate of clustering. 



VI. COMPARISON WITH CLUSTER HALO 
ABUNDANCES 

Since the Liouville cosmology appears not to be ruled 
out by the growth rate data, it is important to define 
observational criteria that will enable us to distinguish 
between it and the concordance ACDM cosmology. An 
obvious choice, that has been extensively used so far, is 
to extract the theoretical predictions of the models for 
the cluster-size halo redshift distributions and to con- 
front them with the data. This will allow distinguish- 
ing the Q-cosmology from the ACDM model. Recently, 
the halo abundances predicted by a large variety of DE 
models have been compared with those corresponding to 
the ACDM model 64]. As a result, it is suggested that 
many DE models (including some of modified gravity) 
are clearly distinguishable from the ACDM cosmology. 

We use the Press and Schechter [41j formalism, which 
determines the fraction of matter that has formed 
bounded structures as a function of redshift. Mathe- 
matical details of such a treatment can be found also 
in [64j; here we only present the basic ideas. The number 
density of halos, n(M,z), with masses within the range 
(M, M + 8M) is given by: 

n(M,z)dM = JL^^f( a )dM , (38) 

where p is the mean mass density. In the original 
Press-Schechter (PSc) formalism, /(<r) = /ps c (cr) = 
y / 2/n(S c /a) exp(— (5^/2cr 2 ), S c is the linearly extrap- 
olated density threshold above which structures col- 
lapse [107], while a 2 (M,z) is the mass variance of the 
smoothed linear density field, extrapolated to redshift 
z at which the halos are identified. It depends on 
the power-spectrum of density perturbations in Fourier 
space, P(k) , for which we use here the CDM form ac- 
cording to [108J . and the values of the baryon density 
parameter, the spectral slope and Hubble constant ac- 
cording to the recent WMAP7 results [6j. Although the 



original Press-Schechter mass-function, /ps c , was shown 
to provide a good first approximation to that provided 
by nuinerical simulations, it was later found to over- 
predict/under-predict the numb er of low/high mass halos 
at the present epoch [l09l IllOJ ] . More recently, a large 
number of works have provided better fitting functions 
of /(<r), some of them based on a phenomenological ap- 
proach. In the present treatment, we adopt the one pro- 
posed by Reed et al. [42j |. 

We remind the reader that it is customary to 
parametrize the mass variance in terms of as, the rms 
mass fluctuations on scales of 8/i _1 Mpc at redshift z = 0. 
In order to compare the mass function predictions of 
the different cosmological models, it is imperative to use 
for each model the corresponding value of cr 8 and 5 C . 
It is well known that for the usual ACDM cosmology 
8 C ~ 1.675, while Weinberg & Kamionkowski [53[ pro- 
vide an accurate fitting formula to estimate 5 C for any DE 
model with a constant equation of state parameter (see 
their Eq.18). Although these conditions are not satisfied 
by our models, one can show [6l|, [6J] that the S c values 
for a large family of dark energy models with a time- 
varying EoS parameter can be well approximated using 
the previously discussed fitting formula, as long as the 
EoS parameter remains constant near the present epoch. 
Since for the current Liouville cosmological model {Q% 
with ft m = 0.28) the effective dark energy EoS parame- 
ter remains close to wde — —0.90, prior to the present 
time, implies that we can use the approximate Weinberg 
& Kamionkowski |53| fitting formula. Doing so we find 
S c ~ 1.668. As a further consistency check we obtain the 
S c value using the same methodology with that of Pace 
et al. (5U; see their Eq. (18)). We find S c ~ 1.676, which 
is in relatively good agreement with that value found by 
the fitting formula of [5J]. Note that for Q m = 0.3 we 
obtain S c ~ 1.674 
are based on the notation of Pace et al. (6l| 

Now, in order to estimate the correct normalization 
of the model's as power spectrum, we use the formula- 
tion presented in Ref. [64[ which scales the observation- 
ally determined as, a value to that of any cosmological 
model. The corresponding value for the Liouville model 
is a s .Q 2 = 0.794 and it is based on ct 8j a = 0.811 (as indi- 
cated also in Table HV) . based on the WMAP7 results [6], 
but also consistent with an average over a variety of re- 
cent measurements (see also the corresponding discussion 
in [64|). For the same model but with fl m = 0.3 we ob- 
tain as, Q 2 = 0.847. 

Given the halo mass function from Eq. (|3"5)) we can 
now derive an observable quantity which is the redshift 
distribution of clusters, Af(z), within some determined 
mass range, say Mi < M//i _1 M < M 2 = 10 16 . This 
can be estimated by integrating over mass the expected 
differential halo mass function, n(M, z): 



Below we utilize those 5 C values which 



M(z 



dV 

dz 



M-2 



n(M,z)dM, 



(39) 



A/, 
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TABLE IV: Numerical results: The 1 st column indicates the cosmological model. The next three columns list the corresponding 
fl m , 0"8 and 6 C values, respectively. The remaining columns present the fractional relative difference of the abundance of halos 
between the Liouville and the ACDM cosmology for two future cluster surveys discussed in the text and in various redshift 
intervals. Error bars are 2<r Poisson uncertainties and are shown only if they are larger than 1CP 2 . 



Model 


* 'm 


C"8 


S c 


(<5A/'/A/'a)orosita 




(<5A/7A/a)spt 












z < 0.3 0.6 < z < 0.9 


z <0.3 


0.6 < z < 0.9 


1.2 < z < 2 


ACDM 


0.28 


0.811 


1.675 


0.00 0.00 


0.00 


0.00 


0.00 


Liouville Q2 


0.28 


0.794 


1.676 


-0.05 -0.27 ±0.01 


-0.06 


-0.15 


-0.37 


ACDM 


0.30 


0.845 


1.675 


0.00 0.00 


0.00 


0.00 


0.00 


Liouville Q2 


0.30 


0.847 


1.674 


0.00 -0.16 ±0.01 


0.00 


-0.02 


-0.26 



where dV/dz is the comoving volume element 

dV 2 dr 

— =n s r 2 (z) — 
dz dz 



(40) 



and Cl s is the solid angle. In order to derive observation- 
ally relevant cluster redshift distributions and therefore 
test the possibility of discriminating between the Liou- 
ville and the ACDM cosmological models, we will use 
the expectations of two realistic future cluster surveys, 
extensively used in recent literature: 



(a) the eROSITA satellite X-ray survey, with a flux limit 
of /iim = 3.3 x 10~ 14 ergs s _1 cm~ 2 , at the energy 
band 0.5 - 5 keV and covering fl s ~ 20000 deg 2 of 
the sky; 

(b) the South Pole Telescope SZ survey, with a limiting 
flux density at i>q — 150 GHz of fu ,iim = 5 mJy and 
a sky coverage of fl s ~ 4000 deg 2 . 

To realize the predictions of the first survey we use the 
relation between halo mass and bolometric X -ray lumi- 
nosity, as a function of redshift, given in Ref. [lll| : 



L(M,z) = 3.087 x 10 



it 



ME{z) 



lO^h-^Mr. 



1.554 



ergs 



(41) 

The limiting halo mass that can be observed at redshift 
z is then found by inserting in the above equation the 
limiting luminosity, given by L — 47rd 2 /i; m Cf,, with dh the 
luminosity distance corresponding to the redshift z and q, 
the band correction, necessary to convert the bolometric 
luminosity of Eq. (jUJ to the 0.5-5 keV band of eROSITA. 
We estimate this correction by assuming a Raymond- 
Smith (1977) plasma model with a metallicity of 0.4Z Q , 
a typical cluster temperature of ~ 4 keV and a Galactic 
absorption column density of uh — 10 21 cm~ 2 . 

The predictions of the second survey can be realized 
using again the relati on be tween limiting flux and halo 
mass obtained in Ref. Ill ill: 



Jvo,}im 



2.592 x 10 8 mJy 

dW) 



M 



W)^h- l M G 



1.876 



E^(z) , 
(42) 
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FIG. 4: The expected cluster redshift distribution of the Li- 
ouville Q2 model and ACDM (dot-dashed curve) models for 
Q m = 0.28 and for the case of two future cluster surveys (up- 
per panels). The lower panels show the corresponding frac- 
tional difference with respect to the reference ACDM model 
(short-dashed curve), while the corresponding curve for the 
case of Q m — 0.3 (see Table ITV| is shown as the long-dashed 
curve. Error bars are 2<r Poisson uncertainties. 



where g?a(-z) = di/(l ± z) 2 is the angular diameter dis- 
tance out to redshift z. 

Below we shall provide predictions for the Liouville 
cosmology only for the Qi model (see Table QJ. For 
the other Liouville models we find that our results re- 
main practically the same (due to the small differences 
in the corresponding Hubble flows). In Fig. [4] (upper 
panels) we present the expected redshift distributions 
above a limiting halo mass, which is M\ = Mri m it — 
max[10 13 /i _1 M0,M/], with Mf corresponding to the 
mass related to the flux-limit at the different redshifts, 
estimated by solving Eq. flU]) and Eq. (g^J) for M, In the 
lower panels we present the fractional difference between 
the Liouville model and ACDM. The error-bars shown 
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correspond to 2<t Poisson uncertainties, which however 
do not include cosmic variance and possible observational 
systematic uncertainties, that would further increase the 
relevant variance. 

The results (see also Table IIV|) indicate that the red- 
shift variation of the differences between the Q cosmology 
and ACDM model is only slightly affected by variations 
in the value of Vt m . For the best- fit case of Q m = 0.28 we 
find that significant model differences should be expected 
for z > 0.6, with the Louville model abundance predic- 
tions being always less than those of the corresponding 
ACDM model. Therefore, we have verified that there 
are observational signatures that can be used to differen- 
tiate the Liouville model from the ACDM and possibly 
from a large class of DE models (see j64[). In Table HVl 
one may see a more compact presentation of our results 
including the relative fractional difference between the 
Liouville model and the ACDM model, in characteristic 
redshift bins and for both future surveys. 



VII. CONCLUSIONS 

In this paper, we have studied the overall dynamics of 
the Liouville cosmological model in which the dark en- 
ergy component is associated with the dilaton. In this 
framework, we first performed a joint likelihood analysis 
in order to put tight constraints on the main cosmolog- 
ical parameters by using the current observational data 
[SNIa, BAOs and H{z)] . We found that the above models 
can accommodate a late-time accelerated expansion. 

Secondly, we performed an additional detailed statisti- 
cal analysis based on the observed growth rate of cluster- 
ing and found that the predicted growth rate of various 
Liouville cosmological models match well the observed 
growth rate. 

The redshift-dependence of the halo abundances, pre- 
dicted by the Liouville model, is significantly different 
from that of the flat ACDM cosmology, at relatively high 
redshifts (z > 0.6). Future cluster of galaxy surveys 
can therefore be used to discriminate the Liouville model 
from other contender DE models. 



under the PCI project A/030322/10 and by the grant 
UV-INV-EPDI1 1-42955 of the University of Valencia. 



Appendix A: Matter density perturbations in dark 
energy cosmologies 

As discussed in the text, the Liouville cosmology, ex- 
plored in this paper, is close to GR in the matter domi- 
nated era (see equation (p} and [66]). Hence, in this ap- 
pendix we would like to remind the reader, for complete- 
ness, of some basic elements of the equations that govern 
the evolution of the mass density contrast, modelled as 
an ideal fluid in the framework of non- interacting DE cos- 
mologies that we employ in this work. Also we assume 
here that for small scales (much smaller than the hori- 
zon) the (effective) dark energy component is expected 
to be smooth and thus we consider perturbations only on 
the matter component of the cosmic fluid. 

In what follows we will adopt a non-relativistic descrip- 
tion, based on a continuity equation, together with the 
Euler and Poisson equations. These equations are: 



|f) V, -(,.111 = . I (AD 



^) (u-V,.lu=-V<I». (A2) 



and 



V^$ = 4irG P + 4irGp DE (3w DE + 1), (A3) 



where (r, t) are the proper coordinates, u is the velocity 
of a fluid element of volume, p is the mass density, pde 
is the DE density, wde is the EoS parameter (given in 
our case by Eq. (|27|l) and $ is the gravitational potential. 
Note that we use here P m = and Pde = wdepde- 

Upon changing variables from proper (r, t) to comoving 
(x, t) ones, the fluid velocity becomes 



u = ax + ax = dx + v(x, t), 



(A4) 
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while the corresponding differential operators take the 
following form: 



and 



\7 X = V = aV r , (A5) 



£ = (£) //xV ,A(il 



where x = r/a and v(x, t) is the peculiar velocity with 
respect to the general expansion. Note that the mass 
density is written as 



P = Pm{t)[l + <5 m (x,£)] 



(A7) 
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In this context, using the standard evolution equation for 
the energy density of ordinary dust, 

p m + 3Hp m = 



as well as Eas.(|A5p. (|A6I) and neglecting second-order 
terms (S m <C 1 and d < u) we rewrite Eqs. dAll ). (|A2[) 
and (|A3|) as: 



dv TT V<1> 

ax + — - + Hv = , 

at a 



V- v = -a- 



dS m 
dt 



(A8) 



(A9) 



— V 2 $ = 4irG Pm {l + 6 m ) + 4irG(3w DE + l)pDE- (A10) 

Following the notation of |88| . we can write the gravita- 
tional potential as follows 

2 2 

$ = 0(x, t) + -TrGp m a 2 x 2 + -nG(3w DE + l)pDEa 2 x 2 . 

(AH) 
Thus, utilizing the Friedmann equation and Eqs. (IAll[) . 
we arrive, after some algebra, at 



and 



dv V0 

— +Hv = , 

at a 



V <j) = AnGa p m 8 n 



(A12) 



(A13) 



Finally, by taking the divergence of Eq. (|A12|) and us- 
ing Eqs. (|A9j) and (|A13j) . we obtain the time evolution 
equation for the matter fluctuation field 



2H5 m - 4:TTGp m 5„ 



(A14) 



In this context, 5 m (a) oc D(a), where D(a) is the lin- 
ear growing mode (usually scaled to unity at the present 
time) . Note that for those cosmological models which ad- 
here to general relativity (GR), the quantity G reduces 
to the usual Newton's gravitational constant Gn^ while 
in the case of modified gravity models (see [72 
we have G = G eS = Y(a)G N (with Y{a)^l). 



A well known approximate solution to Eq. (|31[) is found 
by [82j , where a growth index 7 was used to parameterize 
the linear growing mode for any type of DE with a time 
varying equation of state. Specifically, the approximated 
growth factor was defined through 



D{a) = cxp 



a n£(«) 



dv, 



(A15) 



where now we encapsulate any modification to the G in 
the term fl m (a) = |^f = n E r £ a s ■ This completes our 
discussion on the growth of perturbations due to ordinary 
matter components of the Liouville cosmology explored 
in this work. 
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